Pairing symmetry in BiS2— based superconductors 
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The possible pairing symmetries for BiS2— based superconductors is investigated by using a min- 
imal two-orbital model with onsite and nearest-neighbor intraorbital attractions Vo and Vi, re- 
spectively. By using the mean-field approximation and solving the self-consistent equations, the 
phase diagram of the pairing symmetry is obtained. It is shown that the model allows three 
possible pairing symmetries, depending on the values of Vo and Vi: the isotopic s— wave pairing 

[Ak ~ As], the anisotropic s— wave pairing [Ak = As H ^{coskx -\- cosfc^)] and the d— wave 

pairing [Ak = -2^(cosfci: — cosfcy)]. Furthermore the density of states for these pairing symmetries 
exhibit different behaviors which can be used to distinguish them. 

PACS numbers: 74.70.-b, 74.20.Rp, 74.25.-q 



Introduction. — The recently discovered family of 
BiS2— based superconductors has attracted much atten- 
tion due to its similarity with the cuprates and iron 
pnictides. It displays a layered structure where super- 
conductivity is believed to occur within the BiS2 plane, 
similar to the CuO and FeAs planes in the cuprates 
and iron pnictides, respectively. Superconductivity with 
Tc = 4.5K was first reported in Bi404S3 [l[. Later it was 
found that ReOi-:,¥^B\B2 [Re =La Nd, Ce and Pr) 
can also exhibit superconductivity [2|-|5[ with the highest 
Tc = 10.6K reported in LaOo.sFo.sBiSz 0- These find- 
ings suggest that the BiS2— based superconductors can 
also have relatively high transition temperature and it is 
of great importance to understand the superconducting 
(SC) pairing mechanism and symmetry in this kind of 
materials, since studying these may help to unravel the 
mystery of the pairing mechanism in high-temperature 
superconductors . 

The band structure of this kind of materials has been 
calculated by first principles calculation [6|-|8[, where the 
energy bands close to the Fermi level can be reproduced 
by a simplified two-orbital model [6|. It was shown that 
around x w 0.5, the Fermi surface topology changes and 
the good nesting of the Fermi surface in this case may 
be the cause of the high Tc- Meanwhile the SC symme- 
try is predicted to be s— wave with a constant gap sign 
if the electron-phonon coupling is important, whereas a 
sign-reversing s— wave gap can be obtained if the spin 
fluctuation plays the main role in the Cooper pairing. In 
addition, other pairing symmetries have also been pro- 
posed dunl- 
in this paper, we study the possible pairing symmetries 
in a minimal two-orbital model with onsite and nearest- 
neighbor (NN) intraorbital attractive interactions Vq and 
Vi , respectively. By using the mean-field approximation, 
the phase diagram of the pairing symmetry is obtained. 
We found that in this model there exist three possible 
pairing symmetries, depending on the values of Vb and 
Vi. The first one is the isotopic s— wave pairing [Ak = 
As]. The second one is the anisotropic s— wave pairing 



[Ak = A 

d— wave pairing [Ak 



■^^(cosfc^: -I- cosky)] and the last one is the 



Ai(cosfca: 



- cos fc„ )1 . Furthermore 



we propose that the density of states (DOS) for these 
pairing symmetries exhibit different behaviors which can 
be used to distinguish them. 

Method. — We begin with the minimal two-orbital 
model for BiS2— based superconductors, the mean-field 
Hamiltonian can be written as 



H = ^Z^A^kV'k, 



vi = ( 


c]clt'42t'C-kU'C-k2i), 
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-2ti(cos kx + COS ky) — 2^2 cos(fc^ — ky) 
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-2t3COS(fcj; + ky) 
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-2t4[cos(2A;j; + ky) + cos(fc^ -I- 2ky)] 
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-2t^[cos{2k^, - ky) + cos(fc^ - 2ky)l 


esk = - 
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-xyk -^ 


-2t6(cos kx — cos ky) 
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-2t7(cos2fc^ -cos2fey) 
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-4t8(cos 2kx cos ky — cos k^ cos 2ky). 



(1) 

Here Cy^^^ creates a spin t electron with momentum k 
and in orbital 1. /i is the chemical potential, ii • • -tg are 
the hopping integrals and we consider only spin singlet 
intraorbital pairing up to the NN sites, thus the pairing 
order parameters can be expressed as 



A^k = 



A:, 



S/S, 



As^ H (cos krj., + cos ky) 

. Ad 



2 
- (cos kx 



cos k,. 
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with /3 ~ 1,2 being the orbital index and 
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ifcy)(C-kWCk/3t), (3) 



are the isotropic s— , extended s~ and d— wave compo- 
nents, respectively. N is the number of the lattice sites 
and the doping level x is determined through 



N ^^ 

k/3cr 



--k/3cr 



Ck^cr}- 



(4) 



Eq. ([T]) can be solved as follows: First start with a 
set of random A^/j. A^sp, Adp and /z, the Hamiltonian 
is numerically diagonalized. Then the set of pairing 
order parameters and doping level are calculated by 
using Eqs. ([3|) and ([4]) for the next iteration step 
and /i is adjusted according to the desired doping 
level (xd)- The above procedure is repeated until the 
absolute error of the order parameters between two 
consecutive steps is less than 10~^ and |x — x^l < 10"''. 
Thus by varying the values of Vq and Vi, the phase 
diagram of the pairing symmetry can be obtained. In 
the following, the parameters are chosen as ii • • • is = 
0.167, -0.88, -0.094, -0.014, -0.069, -0.107, 0.028, -0.02, 
respectively. The doping level Xd and temperature T 
are fixed at 0.55 and 0, respectively. We vary Vq (Vi) 
from —0.5 to to get the phase diagram of the pairing 
symmetry. 

Results. — The Fermi surface of the two-orbital model 
for the doping level x = 0.55 is shown in Fig. [T]where the 
small pockets around (±7r, 0) and (0, ±7r) emerge when 
X > 0.515. The self-consistently solved order parameters 
satisfy 



Asp = A,, 
Ad/3 = Ad, 



(5) 



for /? — 1,2, thus the subscript f3 is omitted in the fol- 
lowing. In this case, the intraorbital pairing leads to the 
intraband pairing after a unitary transformation and the 
pairing order parameters in the band representation can 
still be written as 



Ak = As H ;^(cosfc2; -f cosfcj,) 

+ —-{coskx — cosky). 



(6) 



The phase diagram of the pairing symmetry as a func- 
tion of Vq and Vi is shown in Fig. [31 We can see that 




FIG. 1; The Fermi surface of the two-orbital model for the 
doping level x = 0.55. 



the onsite and NN pairings compete with each other. 
When I Vol < \Vi\ (the area filled with blue triangles), 
NN pairing wins over the onsite one and the pairing sym- 
metry is d— wave. On the contrary, when |Vb| > |Vi|, 
the onsite pairing is dominant (see the area filled with 
black squares) and the symmetry in this case is isotropic 
s— wave. Interestingly, in a small parameter range shown 
as the area filled with red circles, the onsite and NN pair- 
ings coexist with each other and the symmetry is predom- 
inantly isotropic s— wave with a small extended s— wave 
component, thus we call it anisotropic s— wave pairing. 
Furthermore, for |Vb| < 0.3 and \Vi\ < 0.2, no SC pairing 
can exist at all. The pairing phase diagram can be under- 
stood as follows: the NN attractive interaction Vi favors 
d— wave pairing while the onsite attractive interaction Vq 
leads to isotropic s— wave symmetry. These two pairing 
symmetries exclude each other, therefore there is no co- 
existing region of them. Only when Vq is strong enough 
and \Vi\ decreases, can the extended s— wave pairing be 
induced between the NN sites, thus it must be accompa- 
nied by an isotropic s— wave component and cannot exist 
alone. 

Then we investigate the possible experimental signa- 
ture of different pairing symmetries through the DOS, 
which can be measured by scanning tunneling microscopy 
(STM). The DOS is expressed as 



p{uj) 



-— ^/77l((Ck/3<T|ct 
k/3cr 



'k/3cr//l^+*+- 



(7) 



Here Im{{. . ■))uj+io+ stands for the imaginary part of the 
retarded Green's function. In Fig. |3]we plot our calcu- 
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FIG. 2: (color online) The phase diagram of the pairing sym- 
metry as a function of Vo and Vi. In the area filled with 
black squares, the isotropic s— wave pairing dominates. In 
the area filled with blue triangles, the d— wave component is 
dominant, whereas in the area filled with red circles. As with 
a small amount of A^s wins over, which we denote as the 
anisotropic s— wave pairing. 



lated DOS. For the isotropic s— wave pairing (the black 
dotted Une), a full gap develops in the SC state and a 
single pair of resonance peaks appears at w ~ ±1. For 
the anisotropic s— wave case, we choose A^s ~ O.SAg and 
the DOS (the blue solid line) is similar to the isotropic 
s— wave case. However, around uj ~ ±1, the resonance 
peak clearly splits into two, leading to a two-gap struc- 
ture. Contrastingly, for the d— wave pairing, the DOS at 
w « is finite and shows a linear dispersion, indicating 
the existence of nodes. Thus the different behaviors of 
the DOS for the three pairing symmetries can be used to 
distinguish them by STM experiment. 

Summary. — In summary, we have studied the possi- 
ble pairing symmetries in a minimal two-orbital model 
for BiS2— based superconductors, with onsite and NN in- 
traorbital attractive interactions Vq and Vi , respectively. 
By using the mean-field approximation, the phase dia- 
gram of the pairing symmetry is obtained. We found that 
in this model there exist three possible pairing symme- 
tries, depending on the values of Vq and Vi. The first 
one is the isotopic s— wave pairing [Ak = A^]. The 



second one is the anisotropic s— wave pairing [Ak = 
As -I- ^j^{coskx -f cosfcj,)] and the last one is the d— wave 
pairing [Ak — -^-{coskx — cosfcj)]. Furthermore we 
propose that the DOS for these pairing symmetries ex- 
hibit different behaviors which can be used to distinguish 
them. 
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FIG. 3: (color online) The DOS as a function of the reduced 
energy to, for the isotropic s— (black dot), anisotropic s— 
(blue solid) and d— wave (red dash) symmetries, respectively. 
For the isotropic s— and anisotropic s— wave cases, cu = uj/As 
while for the d— wave case, w = cj/A^. 
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